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Abstract 

A mesoscopic model of a diblock copolymer is used to study the stability of 
a uniform lamellar phase under a reciprocating shear flow. Approximate viscos- 
ity contrast between the microphases is allowed through a linear dependence of 
the (Newtonian) shear viscosity on monomer composition. We first show that 
viscosity contrast does not affect the composition of the base lamellar phase in 
an unbounded geometry, and that it only couples weakly to long wavelength 
perturbations. A perturbative analysis is then presented to address the sta- 
bility of uniform lamellar structures under long wavelength perturbations by 
self-consistently solving for the composition and velocity fields of the perturba- 
tions. Stability boundaries are obtained as functions of the physical parameters 
of the polymer, the parameters of the flow and the initial orientation of the 
lamellae. We find that all orientations are linearly stable within specific ranges 
of parameters, but that the perpendicular orientation is generally stable within 
a larger range than the parallel orientation. Secondary instabilities are both 
of the Eckhaus type (longitudinal) and zig-zag type (transverse). The former 
is not expected to lead to re-orientation of the lamella, whereas in the second 
case the critical wavenumber is typically found to be along the perpendicular 
orientation. 
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1 Introduction 



Diblock copolymers are macromolecules comprising two chemically distinct and mutu- 
ally incompatible segments (monomers) that are covalently bonded. The equilibrium 
properties are determined by the degree of polymerization A^, (i.e., the length of 
the polymer chain), the volume fraction of one of the monomers /, and the Flory- 
Huggins interaction parameter between the distinct segments x H) 0- While the 
degree of polymerization and the monomer volume fraction are determined by the 
processing conditions, the value of the parameter x is entirely determined by the 
choice of monomers and temperature. 

Above the order-disorder transition temperature Tqdti the equilibrium phase is 
disordered and the monomer concentration uniform. In mean field theory, the order- 
disorder transition takes place at — 10 5 where the Flory-Huggins parameter and 
temperature T are related through x = a/T+P, where a > and (3 are two constants 
m. Below ToDT, equilibrium structures of a wide variety of symmetries have been 
predicted and experimentally observed 0. Around / = 0.5 (symmetric mixture), a 
so called lamellar phase is observed, in which nanometer sized layers of A and B rich 
regions alternate in space. When the copolymer is quenched from a high temper- 
ature to a temperature T < Tqdt, a transient polycrystalline configuration results 
comprising many lamellar domains (or grains), with lamellar normals of arbitrary 
orientations. In practice, full development of the equilibrium state requires very long 
annealing times until substantial long ranged order at the scale of the system size 
can be achieved. However, the underlying ordering mechanisms and associated rates 
that contribute to the large scale reorientation of the grains are essentially unknown 
at present. Not only partial ordering is detrimental for some applications, but it can 
also lead to aging of the material, as well as to potentially anomalous response to 
applied stresses. We present below our analysis of one of the possible mechanisms 
contributing to the kinetics of large scale reorientation of lamellar domains, especially 
in connection with the use of reciprocating shears to accelerate grain coarsening. 

Imposing a reciprocating shear is one of the methods currently in use to achieve 
long-ranged order of block copolymer microstructures. Early work on the response 
of block copolymer blends to shears ^ ^ aimed at establishing the dependence 
of Tqdt on the shear rate, but the experiments also revealed that the shear helped 
select specific lamellar orientations. These observations have subsequently led to a 
large number of groups attempting to quantify the type and degree of ordering that 
can be achieved by reciprocating shears. Early experiments by Koppi et al. 
involved the copolymer poly- (ethylene-propylene) - poly-(ethylethylene) (PEP-PEE). 
Upon lowering the temperature below Tqdt, they observed a transition to the so- 
called parallel lamellae at moderate shear rates, but also an unexpected transition 
to perpendicular lamellae at high frequencies (in parallel alignment, the layers are 
normal to the shear gradient direction, whereas in perpendicular alignment the layers 
are normal to the vorticity direction; see a schematic representation in Fig. |l|) . It was 
also possible to induce the transition from parallel to perpendicular by increasing the 
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shear frequency, but this transformation was not reversible. Koppi et aL interpreted 
the high frequency behavior as shear disordering of the original configuration, followed 
by formation of the perpendicular orientation. Instead, the parallel orientation was 
argued to result from defect mediated growth. 

This phenomenology is qualitatively consistent with the most complete theoretical 
analysis to date due to Fredrickson 0, although one must note that his results were 
obtained for steady shears instead. He used the same model equations that we use 
below in our study, but explicitly allowed for thermal fluctuations near Tqdti as he 
was primarily interested in modeling orientation selection at Tqdt through anisotropic 
fluctuation suppression by the shear flow. He approximated the effect of the flow by 
introducing a modified order parameter mobility that depended on the integrated 
flow over the volume of the sample. Given the inverse characteristic decay time of 
concentration fluctuations, 7*, he showed that for shear rates 7 < 7* the parallel 
orientation is preferred. In the opposite limit, the first transition upon lowering the 
temperature leads to a perpendicular orientation. Further decrease in temperature led 
to transition back to a parallel structure. The range of existence of the perpendicular 
orientation was argued to decrease with increasing viscosity contrast between the two 
microphases. 

However, the phenomenology just described is exactly reversed for a poly-(styrene) 
- poly-(isoprene) (PS-PI) copolymer. Here the perpendicular orientation is observed 
at low frequencies, and the parallel orientation at large frequencies P, There is 
disagreement in the results at low frequencies, as Wiesner and coworkers find a parallel 
orientation at low frequency ||10|, |ll| in the same range investigated by others. The 
experimental results have been summarized in |T^: for frequencies < a parallel 



orientation has been found in two out of four studies. At intermediate frequencies 
ujd < uj < oj'^ the preferred orientation is perpendicular, whereas for uj > uj'^ the 
observed orientation is parallel. The frequency ujd is a characteristic inverse time 
of local domain deformation, and uj'^ is a frequency above which chain relaxation 
dynamics dominates the storage modulus G'{uj). 

A different line of theoretical investigation has shifted the focus of study away 
from fluctuations at Tqdt and into secondary instabilities of a well developed lamel- 
lar pattern [0, 0, Kodama and Doi used a cell dynamical model to study 



possible instabilities of a lamellar pattern upon shearing. The numerical results ob- 
tained motivated in turn an analytic stability analysis that, now in the absence of 



flow, addressed lamellar stability against a change in wavelength. Shiwa |T^ later 
investigated the similarity between the amplitude or envelope equations describing 
slow modulations of a lamellar structure and the same equations describing roll dis- 
tortions in Rayleigh-Benard convection. In the limit of vanishing shear amplitude, 
he established the equivalence of the stability diagram of both systems, and hence 
inferred the range of stability of a lamellar phase against Eckhaus (longitudinal) and 
zig-zag (transverse) instabilities. We were later able to obtain lamellar solutions un- 



der uniform steady and oscillatory shears of finite amplitude |T^, and to study their 



stability. Stability boundaries of the base lamellar phase against both melting and 
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long wavelength perturbations were obtained for transverse and parallel orientations. 
However, we did not address the perpendicular orientation, nor the effect of viscosity 
contrast on the stability diagram. The consideration of both a full three dimensional 
geometry and viscosity contrast are the subject matters of this paper. 

We present a full three dimensional stability analysis of a uniform lamellar phase 
under a reciprocating shear flow with an assumed linear dependence of the shear 
viscosity on local monomer composition. We self-consistently determine the velocity 
field and monomer composition and study the growth or decay of long wavelength 
perturbations of a base lamellar phase. The stability analysis leads to a Floquet prob- 
lem for the perturbation amplitudes, which we solve numerically. First, we find that 
viscosity contrast has a negligible effect on the stability boundaries of the lamellar 
phase. Second, from the actual computation of these boundaries we find that, in 
general, the region of stability is larger for the perpendicular rather than the parallel 
orientation. Third, the marginal mode for the transverse instability branch is typ- 
ically along the perpendicular orientation, so that the initial decay of both parallel 
and transverse uniform states would lead to the appearance of a perpendicular com- 
ponent. We finally show that the orientation of the marginal mode of instability can 
be generally understood from geometrical considerations. At a zig-zag boundary, for 
example, the marginal mode tends to be oriented along the direction that causes the 
largest decrease in lamellar wavelength by the shear. 



2 Mesoscopic model equations and lamellar phases 

At a mesoscopic level, a block copolymer melt is described by an order parameter ipi^^) 
which represents the local density difference between the two monomers constituting 



the copolymer. The corresponding free energy was derived by Leibler [16| in the 



weak segregation limit (close to Tqdt), and later extended by Ohta and Kawasaki 



17 | to the strong segregation regime. If the temporal evolution of ip occurs through 
advection by a flow field as well as through local dissipation driven by free energy 
reduction, ip obeys a time dependent Ginzburg-Landau equation |^ , 

^ + (v • V)ip = MV^i-rijj + uifj^ - KV^i)) - MBip. (1) 

with V the velocity, B a parameter that characterizes the amplitude of the long 
ranged interactions arising from the covalent bond connecting the two sub-chains 



171 , and M a phenomenological mobility or Onsager coefficient The remaining 
parameters appearing in the equation can be related to physical properties of the 
polymer chains as follows |jl3|: K = &^/3, where b is Kuhn's statistical length of the 
chain, r = {2xN - 7.2)/A^, and B = UA/N%'^. 

Under the assumption of Newtonian behavior, the equation governing the velocity 
field is an extended Navier-Stokes equation for an incompressible fluid. 



9v 



dt 



PoT^+Po(v-V)v = V- Vv + VvM -Vp+— VV", V • v = 0, (2) 



5^. 



5ij 
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where po is the (constant) density, 77 the shear viscosity of the fluid which may depend 
on ip, p is the pressure, and appropriate boundary conditions for both ip and v must 
be introduced. The last term on the right-hand side of Eq. (|^) is required to ensure 
that there cannot be free energy reduction by pure advection of ip fl^. This term 
leads to the creation of rotational flow by curved lamellae that is directed toward 
their local center of curvature. 



We first introduce dimensionless variables following ref. |T5[. Lengths are scaled 
by jK/r which is proportional to the lamellar wavelength, time by K/Mr"^^ the 



characteristic monomer diffusion time, and order parameter ip by \Jr /u. As described 
in that reference, only one dimensionless group remains in Eq. (|1]), B' = BK/r^. In 
dimensionless units (we also omit the prime in B\ it is assumed to be a dimensionless 
coefficient in what follows), Eq. (|I]) reads, 

^ + (v- V)V^ = V'(-V^ + ^'- VV)-5V'- (3) 
at 

The mean field order-disorder transition occurs at i?c = 1/4 with a critical wavenum- 



ber gc = \/l/2. 

Equation (j^) can be further simplified by noting that the order parameter dif- 
fusivity (proportional to of Mr) is much smaller than the kinematic viscosity ///po, 
and that under typical experimental conditions poud/r] ^ 1 as well, where u is the 
angular frequency of the oscillatory shear, and d the thickness of the block copolymer 
layer. We therefore adopt a creeping flow approximation according to which the flow 
field instantaneously relaxes to that determined by the instantaneous configuration 
of the order parameter ip (and the no-slip boundary conditions). We further neglect 
in our present study the term ^V^p in the right hand side of Eq. (0). Therefore 
the only dependence of this equation on the order parameter ip enters though the 
shear viscosity 77 = Ti{ip), as we focus on the effect of viscosity contrast between the 
microphases on the stability of a lamellar structure. 

An important restriction of our calculations is that the fluid remains Newtonian. 
This is in line with previous calculations on this system, but it is an approximation 
that needs to be removed in future work, as already noted in the introduction. We 
consider here a specific linear dependence of the shear viscosity on ip, 

V = Vo + ^1^, (4) 

where rji does not have to be small compared to 779. Within the approximation stated, 
and in dimensionless variables, Eq. (||) reduces to, 

= Vp + V [(1 + r/i^) (Vv + Vv^)l , (5) 



where the dimensionless viscosity correction ri[ = ^Jr /u{rii/riQ) has been introduced 
(and the prime removed). 
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The physical system under consideration is a layer of block copolymer, unbounded 
in the Xi and X2 directions, and being uniformly sheared along the X3 direction (Fig. 
p. The layer is confined between the stationary = plane, and the plane x^ = d 
which is uniformly displaced parallel to itself with a velocity Vpiane = 'jdu cos(ci;t)xi, 
where xi is the unit vector in the xi direction. 7 is the dimensionless strain ampli- 
tude and uj is the angular frequency of the shear. In what follows, this velocity is 
also expressed in the dimensionless variables given above. In particular, the dimen- 
sionless wall velocity can be written as Vp^^^j^g = 7' cos{u't') so that with the scalings 
introduced earlier, 7' = \/2{ujTD){qcd)'y, with td = K/Mr'^ the (diffusion) time scale 
introduced earlier. Since for a realistic system q^d ^ 1, the relatively small values 
of 7' required for instability (see below) imphes that our study addresses in practice 
shear frequencies much smaller than the inverse diffusion time. Again, in what follows 
we restrict ourselves to dimensionless variables and drop the primes. 



We first summarize the results of ref. [jT3[ in which the viscosity was assumed 
uniform. We focused there on the weak segregation limit e = {Be — B)/2Bc ^ 1, in 
which the solution for the monomer composition can be obtained perturbatively in e, 

^{r) = 2A{t) cos(q ■ r) + Ai{t) cos(3q ■ r) + . . . (6) 

where r = xiXi -|- X2X2 + 0:3(7 sin (co't)xi + X3) so that it can be thought of as having 
components in a non orthogonal basis set which follows the imposed shear, and q = 
(51,52,93) is the wavevector in the corresponding reciprocal space basis set {gi = 
Xi— 7 sin(a;t)x3, g2 = X25 gs = X3}. Note that in this new coordinate system, perfectly 
ordered configurations are stationary. Three orientations relative to the shear can be 
defined as follows: 53 7^ 0, gi = 53 = is a purely parallel orientation, 52 7^ 0, gi = 53 = 
is a perpendicular orientation, and gi 7^ 0, 52 = 93 = is a transverse orientation. 

In the absence of viscosity contrast, the velocity field is given by, 

v" = 'juj cos(co't) X3 Xi. (7) 

Furthermore, by substituting Eq. (|^) into Eq. (^, the lowest order solution (0(e^/^)) 
is given by 

V^(r) = 2A(t)cos(q-r) (8) 
where the amplitude A{t) satisfies the equation ||T5[| , 

dA 

— = a[q\t)]A-3q\t)A\ (9) 

with g^(t) = ql + [7sin(a;t)gi — 53]^ + ^2 and o"(g^) = q^ — q^ — B. This equation 
can be integrated to give the marginal stability boundaries, and the function A{t) 



itself IjT^. From this analysis, a critical strain amplitude 7c, which depends on q, 
was identified such that for 7 < 7c the uniform lamellar structure oscillates with the 
imposed shear, but at 7 > 7c A{t) decays to zero; i.e., the lamellar structure melts, 
to use the terminology used in experiments. 



6 



The stability of the base lamellar pattern was then addressed in two spatial di- 
mensions by a Floquet analysis. The study was restricted to the (xi,X3) plane (i.e., 
52 = 0), and therefore restricted to transverse (gs = 0) and parallel (gi = 0) lamel- 
lae only. Briefly, for sufficiently small 7 all orientations were seen to retain a range 
of stability in wavenumber, range that becomes narrower along the diagonal in the 
(gi, ga) plane. The range of stability is further reduced with increasing frequency. At 
moderate shear amplitudes (e.g., 7 = 0.4), only fully parallel or transverse lamellae 
remain stable, and there is also a weak frequency dependence. We extend below these 
results to three spatial dimensions and to a fluid with a nonuniform shear viscosity, 
function of the local monomer composition ip. Our aim is to incorporate into the 
stability analysis a different effective rheology of the two microphases. 



3 Stability of a lamellar phase under shear 

We obtain in this section the flow in the melt that arises from the nonuniform shear 
viscosity, and the resulting corrections to the base lamellar solution. We then per- 
form a self-consistent stability analysis of the lamellar order parameter and flow field 
against long wavelength perturbations. 

There is now ample evidence that viscoelastic contrast between the microphases 
affects orientation selection. We wish to investigate here whether such a contrast 
substantially affects the stability of a uniform lamellar phase to shear. This is an ini- 
tial step in attempting to understand the experimental phenomenology. For example, 
and as briefly discussed in Section 0, the qualitative response to oscillatory shears 
in a system such as poly- (ethylene-propylene) - poly-(ethylethylene) (PEP-PEE) is 
qualitatively very different than that of, say, poly-(styrene) - poly-(isoprene) (PS-PI). 
The former prefers the parallel orientation at low frequencies and the perpendicular 
orientation at higher frequencies, whereas the behavior is essentially reversed for the 
latter. As emphasized by Fredrickson and Bates 0, the microphases of the PEP-PEE 
system are well matched mechanically. However, PS is largely unentangled whereas 
PI is entangled. Therefore a large contrast in the relaxation times of the blocks is 
anticipated. Also recent experimental evidence by Winey and coworkers suggests 
that the response of a parallel conflguration to shear can be qualitatively described by 
a three region model. One central region in the vicinity of the covalent bond between 
the A and B chains is relatively stiff, and responds elastically under shear. This re- 
gion is surrounded by two other regions with a largely viscous response to the shear 
as the chains are elongated. Before addressing the more general case of viscoelastic 
contrast between the microphases, we study in this paper a fairly simplistic situation 
in which the blend remains a Newtonian fluid, but one in which the shear viscosity 
depends explicitly on the local monomer composition rj = and hence the actual 
flow fleld inside the melt depends on the order parameter configuration. We follow 
a related study by Fredrickson ^ who incorporated viscosity contrast by computing 
an averaged, effective shear rate for a lamellar configuration, which then renormal- 
ized the mobility coefficient in the time dependent Ginzburg-Landau equation for the 
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order parameter. We do consider here, however, the full calculation of the flow field 
self-consistently, although for practical convenience we restrict our analysis to the 
linear variation given in Eq. (^. In this case, we can solve for the base flow field 
exactly, fact that considerably simplifies the stability analysis. 

We begin by writing Eq. in components, 

dip - + ^xt\)){diVj + d^Vi)\ = 0, (10) 

and consider periodic boundary conditions along the ici and X2 directions. To facili- 
tate the computation, we replace no shp boundary conditions in the X3 direction by 
sheared periodic boundary conditions (i.e., periodic boundary conditions in a frame of 
reference rigidly attached to the moving plates. See, e.g., ref. |jl^ for further details). 
We then introduce the decomposition 

v = v° + u. (11) 

The governing equation for the velocity u is (with '?/'(r) = 2Acos(q ■ r)) 

- dipi + d,[{l + vi^){div'^ + d,v^)] + djlil + Vi^){diUj + d,u,)] = (12) 

In order to solve this equation, we introduce a new set of coordinates {x',y',z') 
such that z' = q (the z' direction is parallel to the time dependent wavevector of 
the lamellar phase), and y' is perpendicular to the plane (g^, 0, g^)- Then Uy' and u^' 
vanish in the lamellar phase, and the equation for Ux', after some algebra, is given by, 

dz' [(1 + 2riiA{t) cos qz')dz'Ux'] = '-)r]iA{t)q(3 sin qz' , (13) 

with /3 = 2[-{ql - qlf - q^{ql + ql)\/qnx and = {{q^iql - q^ - qD^^qxqyq^, 
~1x{(ll + Qy ~ ^z))!; where the wavevector q is a function of time as the orientation 
of the base lamellar structure adiabatically follows the imposed shear flow. Equation 
([I3D can now be integrated twice to obtain Ux'- It is important to note that for a 
purely perpendicular base state (3 = 0, and hence the flow correction u = in this 
case. This is in agreement with the results of Fredrickson [0. 

A general property of the solution for the flow perturbation u that follows from 
the linear form of the viscosity contrast is that the velocity is always perpendicular to 
q, i.e., parallel to the lamellar planes which are the planes of constant ip. Therefore 
the advection term u ■ Vip = for a lamellar phase, and hence the base fleld ip is 
unaffected by the flow. Note that this result is general as long as the order parameter is 
a function of one spatial direction only. In the more general cases, only the functional 
dependence in cos(z') and sin(z') in Eq. (pTsD would be different. In short, Eq. (§) 
remains the solution for the order parameter when the viscosity contrast is given 
by Eq. (|^). The base flow fleld can be obtained from the solution of Eq. (|13D and 
Eq. m). 



We next address the stability of the solution Eqs. (|^) and (H) against long wave- 
length perturbations at fixed e . We consider solutions of the form 



4j = tpl+tp2 

V = v° + u + w 
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where ipi is the base solution given in Eq. (||), u the solution of Eq. (]T3|), and, 

^2 = 1P2+ e'^i'+'Q" + il)2- e''^'-''^ ' + c.c. (14) 

If the wavevector of the perturbation Q is parallel to the base wavevector q, the 
perturbation is said to be of the Eckhaus type. In the case of Q _L q, we are consider- 
ing a zig-zag perturbation. Due to our choice of shear periodic boundary conditions, 
and unlike more standard analyses of long wavelength instabilities, the perturbation 
wavevector Q is time dependent and periodic with the same periodicity as the shear 
(the same is true for q). 

The stability of periodic solutions of the current model in the absence of flow has 



been extensively studied ||2^, There is a range of wavenumbers q within which 
the lamellar states are stable. This range of stability increases with the distance to 
threshold e. At large wavenumbers, lamellae undergo an Eckhaus instability which 
tends to lowers the value of q by eliminating lamellar layers. On the other hand, 
at small wavenumbers the wavenumber is increased via a zig-zag instability. As we 
show below, the imposed oscillatory shear has important consequences for these two 
instabilities, especially for the zig-zag case. In particular, we show that the way in 
which a zig-zag instability leads to a readjustment of the lamellar wavelength depends 
strongly on the relative orientation between the base lamellae and the shear direction. 

In order to derive an evolution equation for the two amplitudes il)2+ and 'ip2- 
we need the velocity field w. Since the convective term in Eq. (||) is v ■ Vip = 
{u + w)-'V{ipi + ip2), we only need to retain Fourier components of the form e±«q-i"±«Q r_ 
Therefore one only needs the Fourier components e^^*'*''" for u and e^^*'i''"='=**^''" for w. 
The governing equation for the second order velocity w reads, 

-diP2 + dj[{l + r]ii)i){diWj + djWi)] 

+ dj [riiil)2 {diUj + djUi)] + dj [ijiipi (diV^ + djVi)] = 0. 

This equation is solved by transforming again to the (x', y', z') coordinate system. 
The sum of the third and fourth terms is the inhomogeneous term in the equations 
and is of the form of \\[qz')e^^ ^ + c.c. Therefore we write the second order velocity 
as 

w = e*^'''"'w(g2;') + c.c. 
P2 = e**^''''p2(5'-2') + c.c. 

which leads to the equations 

- diP2{qz') + (zQ;. + d,)[l + r/i^i(gz')]x 

[{tQ[ + di)w,{qz') + {iQ'^ + d,)w,{qz')] = h,{qz'). (15) 

We now solve for the functions w(gz') and P2{qz') numerically together with the 
incompressibility condition, 

iQx'Wx' + iQy'Wyi + [iQz' + dz')wz' = 0. (16) 
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Since the velocity w is proportional to 1^2, and given the latter's decomposition 
in its components e^*'' '" and we find it useful to express w as 



Also by considering the Fourier expansion of u, 



w e 



e'^' + c.c. 



c.c. 



we arrive at the (hnearized) system of equations obeyed by the perturbation ip2, 



d 


i^2+ 




di 







H2l{t) 



Hl2{t) 
H22{t) 



^2+ 
^2- 



(17) 



This is the central result of this section. The matrix elements are given by, 



Exx{t) = -1+ ~ ll + B + 6A{t)% + iA{t)w++ ■ q 

Huit) = 3A(t)2/+ + iA(t)w-+-q + m2-(q-Q) 

Hu{t) = 3A{t)H^-iA{t)w+- ■q-iu2-{q+Q) 

H22{t) = -l--t+B + 6A{tfl^-iA{t)w— -q 

with l± = — (q+Q)^. First we note that the matrix elements are periodic in time, with 
the same periodicity as the imposed shear. Therefore the stability analysis constitutes 
a Floquet problem for the amplitudes 1IJ2+ and ^^2- which we solve numerically. Briefiy, 
the function A{t) is evaluated numerically for a given lamellar orientation and shear 
fiow parameters. The fiow correction due to viscosity contrast u in the base state is 
known analytically, but the fiow perturbation w is obtained numerically by solving 
Eqs. ( ]T5D and (|T6|). Once the matrix elements in Eq. (|T8|) have been numerically 
evaluated, the Floquet problem is recast as an eigenvalue problem, which is then 
solved numerically. For fixed e, 7, u and rji the stability boundaries are determined 
as the loci of q at which the eigenvalue function (t(Q) changes from a maximum to a 
saddle point at Q = 0. 

We finally note that the fact that the Floquet problem only involves the amplitudes 
ip2+ and ip2- and not w++, , w ^ and w is a direct consequence of the creeping 
fiow approximation introduced: the velocity field is slaved to the composition field. 



4 Results and discussion 

As illustrated in Fig. |I] and discussed in Section |I], the configuration considered is an 
unbounded layer being sheared along the Xi direction, and with a velocity gradient 
in the X3 direction (and hence the vorticity vector is along the X2 direction). Three 
particular orientations of the lamellae are of special concern in relation with the 
issue of orientation selection under shear: parallel, perpendicular, and transverse, 
as indicated schematically in Fig. |1|. These three special cases are to be discussed 
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first, although we later present the results of the stability analysis for an arbitrary 
orientation of the lamellae. 

We again note that we do not focus on the formation of a particular orientation 
from an initially disordered configuration, but rather on the simpler problem of as- 
certaining the region of stability of a base lamellar state of a given orientation as a 
function of the shear rate. Of course, linearly unstable orientations are not expected 
to be observable in experiments, not even locally during the transient evolution of 
polycrystalline configurations. Our analysis, however, does not address the deter- 
mination of the basin of attraction of each orientation from an initially disordered 
configuration, and hence cannot completely answer the orientation selection problem. 
However, as we will show below, the regions of stability against long wavelength per- 
turbations are quite small, and therefore our results do provide some guidance for the 
orientation selection problem. 

One of the main conclusions from our numerical study is that velocity field correc- 
tions due to the assumed viscosity contrast between the two microphases phases have 
a negligible effects on those stability boundaries that we consider. We believe that 
this is due to two reasons. First, we have shown that viscosity contrast of the assumed 
form has no effect on the base state of perfectly parallel lamellae (c.f. Eq. (p!3| ) and 
the paragraph that follows). Second, both Eckhaus and zigzag instabilities are long 
wavelength instabilities. Consider, for example, a transverse zig-zag perturbation of a 
parallel state. Since the induced flow u alternates in direction on consecutive parallel 
planes, it will distort the zig-zag perturbation, but the effect will be very small for a 
long wavelength perturbation. Hence in what follows we concentrate the discussion 
on the effects that follow from lamellae orientation, keeping in mind that they are 
quite independent of the value of rji considered. 

The stability of a particular lamellar phase (with a base wavevector q) is de- 
termined by the growth or decay of perturbations of wavevector q + Q which is in 
turn given and by the sign of the corresponding eigenvalues of the Floquet problem 
defined by Eq. (p^Tf). Also, the orientation of the marginally unstable mode Q can 
be used as an indicator of the orientation of the emerging structure. Of course, in 
an extended sample in which multiple orientations coexist, possibly with a subset 
of them becoming linearly unstable, it will be their nonlinear competition that will 
determine the asymptotically selected orientation. Nevertheless, we believe that it is 
still useful to catalog the orientation of the marginal mode in the case of the three 
basic orientations of the lamellae: parallel, perpendicular and transverse. 

Before we turn to the numerical results of the Floquet analysis, we also mention 
that we find it useful to interpret them in terms of a simple geometrical construction 
of how the imposed shear affects the lamellar wavelength. If one were to neglect 
any flow corrections due to viscosity contrast, and if the lamellar base state were to 
follow the base shear flow adiabatically, then a lamellar phase with initial wavevector 
q(t = 0) = (o'l, Q'2) in the Cartesian or laboratory frame of reference, evolves 
according to 

q(i) = ill, <i2, -iqi + gs), (18) 
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with gi,g2 and qs constant and 7 = 7sin(ci;t). We will see that in most cases the 
marginal mode at the zig-zag instability is that which leads to the largest increase in 
wavenumber by the shear. 

First consider initially transverse lamellae, i.e., with q along the xi direction. 
From Fig. |l| it is clear that the shear flow will tilt the layers, and in doing so decrease 
their wavelength. The stability diagram for this case has been obtained by numerical 
solution of the Floquet problem for lj = 0.01 and e = 0.04 and is shown in Fig. ^ 
as a function of the shear amplitude 7. The outside bounding curves are the neutral 
stability boundaries so that only within this range a nonlinear solution for ip exists, 
with an amplitude A{t) given by Eq. (^. The bending of the curves toward smaller 
values of g as 7 increases can be qualitatively understood by noting that the oscillatory 
shear leads to a decrease in the lamellar wavelength. Thus at higher values of 7, the 
region of existence of flnite amplitude solutions shifts toward larger wavelengths to 
compensate for a larger reduction in wavelength by the flow. At large enough 7, 
however, nonlinear solutions cease to exist. 

Only nonlinear solutions within the inner region in Fig. |^ are stable against long 
wavelength secondary instabilities. Both Eckhaus and zig-zag type instabilities occur 
as shown in the flgure. The Eckhaus instability is a longitudinal phase instability (Q || 
q) so that through coupling with the amplitude it leads to an amplitude modulation 
in the same direction as the base periodicity. In general, this instability appears in 
the large wavenumber range of the diagram and is qualitatively interpreted as an 
instability that leads to a decrease in the wavenumber by eliminating lamellar layers 
through the amplitude modulation. On the other hand, in the range of small q, a zig- 
zag instability with Q ± q can lead to an increase in q to return the unstable lamella 
to the stable range. In the absence of shear, perturbations to a base state deflned by 
q = gxi along (5x2 or Qic^, are equivalent. Under shear, however, the degeneracy is 
broken, and we flnd from the numerical Floquet analysis that the marginal mode for 
the zig-zag instability is Q^2, i-e., along the perpendicular direction. In summary, 
the large wavenumber instability is of the Eckhaus type and does not lead to lamellar 
reorientation. The small wavenumber instability is of zig-zag type, and leads to a 
growing component along the perpendicular direction. 

We next turn to a geometrical interpretation of these results based on a rigid 
distortion of the lamellar structure as discussed above. The marginal wavevector for 
an unstable zigzag perturbation can be either q + Q = {q,Q,0) which leads to growth 
of a perpendicular component, or q+Q = (g, 0, Q) that leads to a parallel orientation. 
The effective wavenumber of the distortion in each case is given by 

perpendicular zigzag |q+Qp(t) = !(<?, QjT''?)!^ = (1 + 7^)'?^ + 

parallel zigzag |q+Qp(t) = \{q,0,^q + Q)]"^ = {1 + ^^)q^ + 2^qQ + Q^. 

Under oscillatory shear 7 changes sign and both cases lead to a wavenumber increase, 
although given the sign change, it is not possible to tell which mode is more effective at 
increasing the wavenumber following the instability. This is the only case in which this 
geometric interpretation does not lead to a selected marginal mode consistent with 
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the Floquet analysis. We note, however, that if the shear were stationary instead 7 
increases monotonically, and a zig-zag instabihty along the parallel direction would 
provide for a larger increase in wavenumber. 

We now consider a base lamellar phase oriented parallel to the shear. Figure ^ 
shows the stability boundaries computed numerically. As seen in the figure, secondary 
instabilities are again of the Eckhaus type at large g, and of zig-zag type at small q. 
In this case the location of the boundary is nearly independent of 7, i.e., the shear 
flow has no effect on the Eckhaus instability. This can be understood geometrically 
from Fig. since the parallel orientation is unaffected by the shear, and an Eckhaus 
instability of a parallel orientation would also lead to a growing mode along the par- 
allel orientation. The zig-zag boundary at small wavenumbers also remains constant 
and unaffected by the shear. Nevertheless shear flow does introduce a distinction 
between the two possible zig-zag modes along the perpendicular and transverse di- 
rections respectively. The former, given by q = (0, Q, q), will not be advected by the 
shear according to the Eq. (plsl) . The latter is given by q = {Q,0,q), and it has its 
wavenumber advected as q{t) = {q — jQY + ■ Therefore the shear leads to a net 
wavenumber increase only when the instability is along the transverse direction. The 
results of the Floquet analysis are consistent in this case with the orientation that 
produces the largest wavenumber increase. 

We finally discuss the case of initially perpendicular lamellae. The stability bound- 
aries for this case are also shown in Fig. ^. Again, the Eckhaus boundary is unaffected 
by the shear. This can be understood from Fig. |I] since the perpendicular orientation 
does not couple to the flow. An Eckhaus mode is along the perpendicular direction 
as well, and hence remains unaffected by the shear. On the other hand there is a very 
weak dependence of the zig-zag stability boundary on the shear amplitude, and we 
find that the marginal mode is along the transverse direction. This is also consistent 
with the geometric interpretation given earlier since an imposed shear flow will always 
increase the wavenumber for a transverse zig-zag mode (Q, g, 0), but has no effect on 
a parallel zig-zag mode (0, g, Q\. 

transverse zigzag |q+Qp(t) = |(<5, g, 7Q)P = 9^ + (1 + 7^)<5^ 
parallel zigzag |q+Qp(t) = |(0, g, = g^ + 

Our results for the the secondary instabilities of the three particular lamellar 
orientations can be classified according to the orientation of the marginal mode as 
follows: 

Transverse lamellae — > zig-zag to perpendicular orientation 
Parallel lamellae zig-zag to perpendicular orientation 

Perpendicular lamellae zig-zag to transverse orientation 

If one assumes that the orientation of the marginal mode sets a preference for 
the final orientation of the emerging stable lamellar structure, our results suggest 
that the perpendicular alignment is the preferred orientation under the action of the 
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oscillatory shear, independently of viscosity contrast (within the Newtonian viscosity 
model adopted). Even in the case in which an instability of a perpendicular ori- 
entation leads to a transverse state, a shear flow of small amplitude would drive it 
back to the perpendicular orientation, presumably now with a wavenumber inside the 
stable region. We however need to be remain cautious about the generality of this 
conclusion as a complex dynamical evolution could follow the initial decay after a 
zig-zag instability. For example, our earlier numerical results in ref. showed that 
a zig-zag instability may lead to kink-band formation, and yet to another orientation 
change as the bands themselves become unstable. 

We finally present the stability diagram in full three dimensional q space for a 
base lamellar state of arbitrary orientation relative to the shear. In the absence of 
flow, all lamellae with wavenumbers near the marginal value Qc are stable. Hence 
the stability region is a spherical shell in q space. Under shear, any lamella with a 
wavevector that has a significant component along the transverse direction becomes 
unstable, and the region of stability becomes compressed along the qi direction. For 
sufficiently large 7, the stable region adopts a toroidal shape around the q2-<iz plane 
with a small projection in the qi direction that depends on 7. Although we have 
argued earlier that lamellae with a too small wavenumber will become unstable to a 
zig-zag mode and evolve into a new state oriented along the perpendicular direction, 
this by no means implies that a stable parallel orientation cannot exist (although 
perhaps within a narrow range of wavenumbers). 

As an example. Fig. | shows the complete stability boundaries along both q2-qi 
and q^-qi planes for e = 0.04 and 7 = 1. Note how the stable region in the q2-qi 
plane is significantly larger than in the q-^-qi plane. This is also evidenced in Fig. ^ 
by the asymmetry in the toroidal stability region which is wide near the q2 axis (the 
perpendicular direction) but narrow near the q2 direction (the parallel axis). The 
implication of this result is that if an initially disordered state is comprised of many 
domains locally oriented along arbitrary directions, one would expect that a larger 
portion of the sample would remain in the perpendicular state, as it corresponds to 
the one with the largest region of stability in q space. This argument is distinct 
from that given above involving the orientation of the marginal model following a 
secondary instability of a perfectly ordered lamellar structure, but it suggests the 
same dominant orientation. 

In addition to the detailed numerical computation of the stability region in q 
space, the relative sizes of the stability regions for parallel and perpendicular orien- 
tations may be estimated from the geometric interpretation given above. Consider 
first an almost perpendicular state with a wavevector of a magnitude go with a small 
component Ag in the transverse direction. Under the action of the shear, Eq. ([T8| ) 
gives for the time dependent wavenumber 



Therefore the wavenumber increases with the shear as (Ag/go)^- On the other hand. 
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an almost parallel state has a wavenumber 




q{t) = [Aq^ + {jAq + go) 

and therefore it increases only linearly with Ag/gg. These two relations essentially 
yield the dependence of the width of the toroidal region in Fig. ^ along the transverse 
direction, and therefore the fact that the relative extent of the perpendicular region 
of stability is larger than the parallel region. 

Insofar our results could be used to infer the dominant orientation under shear, 
our conclusions would differ from those of Fredrickson from low frequency ex- 
perimental evidence in PEP-PEE, but agree with some low frequency experimental 
evidence in PS-PI. One must caution, however, that as discussed in ref. |]^, it is 
not clear which of the two systems PEP-PEE or PS-PI has a larger value of rji, or 
how close they are to the assumed Newtonian behavior. Therefore we must consider 
our results to be a baseline against which to compare future work involving true 
viscoelastic contrast between the microphases. 

Finally it is interesting to address a possible reason for the qualitative discrepancy 
between the conclusions of our work and those of Fredrickson's concerning the effect 
of viscosity contrast. Although he discussed the case of steady shear, our results are 
confined to low frequencies, and are only weakly dependent on frequency. Therefore 
the discrepancy in conclusions is not likely to arise from the time dependence of our 
solutions. Briefly, he found a transition from a high temperature disordered state to 
parallel lamellae at low shear rates, and to perpendicular lamellar at high shear rates. 
In the latter case, and upon further decrease in temperature, he predicted another 
transition to a parallel state. The location of this second transition line depended on 
the viscosity contrast between the microphases. Instead, we see no discernible effect 
arising from viscosity contrast. A possible qualitative explanation can be given as 
follows: Our calculation is conducted at an externally imposed shear rate D = Vp/d 
in dimensional units (here and in what follows), where Vp is the velocity of the solid 
boundary. Consider a parallel configuration comprised of only two planar layers with 
uniform shear viscosity ?7o + rji and ?7o — rji respectively. The average shear rate of 
this configuration is given by 



d {v' - 0) d {vp - v') 



2 d/2 2 d/2 



where v' is the (unknown) speed at the boundary between the two layers. In this case, 
the average shear rate is independent of the viscosity contrast. On the other hand 
and following Fredrickson, if the boundary conditions involve an imposed shear stress 
a, the conclusion is different. The shear rate for the base state is now D = a/rjQ, the 
same as in the previous case. However, in the creeping flow approximation considered, 
the shear stress is uniform in the fluid, but the shear rate is different in the two layers, 

DiiVo - Vi) = ^2(110 + rji) = Dr]o, 
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where Di (resp. D2) is the shear rate in the layer of viscosity rjo — rji (resp. rjo + rji). 
The spatial average over the configuration is now, 











h -D2 




2 





,1-5 



+ 



5 



D 1+5^ 



where 6 = Tji/rjQ. Therefore the average shear rate across the layer in a parallel con- 
figuration always increases with viscosity contrast. (In both sets of calculations, the 
perpendicular base state is not affected by the shear flow). Fredrickson's results are 
based on the introduction of a renormalized order parameter mobility that is assumed 
to depend on the average shear rate across the layer. Therefore his results for a parallel 
conflguration are affected by the shear, whereas ours are not. As discussed in Section 
the composition fleld of the base parallel conflguration is not modifled by the flow 
correction arising from viscosity contrast. Finally, although the implications of the 
choice of boundary conditions on the dynamical evolution of partially ordered lamel- 
lar conflgurations are difficult to establish, we note that Fredrickson's work focused 
on fiuctuations near a uniform, disordered state, whereas we have considered pertur- 
bations around a (weakly) nonlinear state of a finite amplitude, saturated lamellar 
structure. 
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Figure 1: Schematic representation of the geometry considered including the shear 
direction, and the three different lamellar orientations discussed in the text. 
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Figure 2: Neutral stability curve and secondary instability boundaries for an initially 
transverse lamella with uj = 0.01 and e = 0.04. 
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Figure 3: Neutral stability curve and secondary instability boundaries for initially 
parallel and perpendicular lamellae with lu = 0.01 and e = 0.04. 
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Figure 4: Stability regions in the qi-q2 and qi-qs planes with cu — 0.01, 7 = 1, and 
e = 0.04. 
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Figure 5: Region of stability of a uniform lamella as a function of its orientation. The 
axes correspond to transverse (gi), perpendicular (^2), and parallel (gs) orientations. 
The surface has been determined by numerical solution of the Floquet problem defined 
by Eq. (|T^ with u = 0.01, 7 = 1, and e = 0.04 
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